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Abstract 



Large mixing induced fg transitions can explain the deficit in the flux 

^ \ of the low energy atmospheric whereas the solar neutrino deficit can be 



simultaneously explained through MSW transitions fe — > i^t- A combined 
analysis of both these effects is presented. The large Ve — mixing affects 
the MSW transition between V(, and Vr significantly. As a consequence, a 
large region of parameters ruled out by experiment in the two generation case 
is now allowed. The mass hierarchy as well as the mixing pattern required 
arise naturally in a model for neutrino masses proposed by Zee. 
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Two independent sets of experiments point towards the existence of neutrino oscillations. 
The flux of the solar neutrinos is found to be depleted by varied amounts compared to the 
expectations based on the standard solar model in four different experiments namely, Home- 
stake 0], Kamioka 0, GALLEX |[ and SAGE §]. Similarly, the flux ratio of the low energy 
(< GeV) to z/g of atmospheric origin is seen to be smaller than the theoretically expected 
one in experiments by the Kamioka II and 1MB ^ group. Independent measurements by 
the SoudanQ and the Frejus collaborations are not inconsistent with the above but the 
statistical significance of these results is much less than the corresponding one in cases of 
Kamioka and 1MB results. Theoretically, both these deficits could originate from the same 
common source, namely neutrino oscillations. However the range of the relevant parameters 
((mass)^ difference A and mixing angle 6 ) required by the solar and atmospheric neutrino 
data are quite different. 

The solar neutrino deficit can be explained through neutrino oscillations occurring either 
in vacuum or in the presence of solar matter ||T0|. The former requires A ^ 10^^'^ (^V)^, 
sin^26' ^ 0.75 — 1.0 while the latter becomes important if either A ^ (5 — 10) x 10"^ (eV)^, 
sin^ 29 ^ (0.2—0.9) or A sin^ 29 ^ (3—4) x 10^*^ (eV)^ corresponding respectively to the large 
angle adiabatic and non adiabatic Mikheyev Smirnov Wolfenstein (MSW) solutions [ITT 



In contrast, the deficit in the atmospheric muon neutrinos can be explained by means of 
neutrino oscillations |jl2| when parameters assume typical values A ^ 10^^ — 10^^ (eV)^ and 



sm^29 ^ 0.5. 

It is quite clear that a simple two generation picture used in deriving above restrictions 
is inadequate to simultaneously explain the solar and the atmospheric neutrino deficits. A 
natural possibility is to consider more realistic three generation picture. When all three 
neutrino flavors are involved, the vacuum oscillations can account ||13[ for both the deficits 



*The latest, still preliminary result from the SOUDAN-II detector seems to confirm the results 
of Kamiokande and 1MB, albeit with larger statistical errors. 
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if |Ajj| = \mf — rrijl (i, j=l,2,3) ^ (0.5 — 1.2) x 10~^ (eV)^ and at least two of the mixing 
angles are large. On the other hand if |A| < 10~^ eV^, then the matter effects inside the 
Sun become important and they should be taken into account. There are two independent 
possibihties in this case. 

(z) The oscillations between z/e and generate the solar deficit while the reduction in the 
low energy atmospheric neutrino flux could be due to the i/^ oscillating to Ur- This would 
require a large 623 and IA23I ~ 10~^ — 10^^ (gV)^. In this case, the high energy muon 
neutrinos originating in the atmosphere also undergo oscillations as they travel through the 



Earth. Such oscillations have not been found [|1J,|T^]. This negative evidence combined 
with the positive evidence for oscillations in case of the low energy neutrinos [Q,^ strongly 
constrain the relevant parameters and a very narrow range 2 x 10^'^ < IA23I < 0.4 (eV)^, 
0.4 < sin^ 2^23 ^ 0.7 is consistent with all the available information^ ||17|| . 

{ii) The other possibility is to have a large mixing between z/g and z/^ and relatively large | A21I 
to generate the deficit in the atmospheric flux. The MSW solution for the solar neutrino 
can be obtained in this case if A31 falls in the appropriate range. The upcoming high energy 
muons should oscillate in this well but these oscillations will be affected by the 

presence of matter inside the Earth. The values of parameters which are consistent with 
the observations of the low and high energy atmospheric z/^ flux as well with the laboratory 
search for z/g — z/^ oscillations were determined to be fl^: 4 x 10^^ < IA12I < 0.02 (eV)^, 
0.35 < sin2 20i2 < 0.7 . 

In case (i), the solar and the atmospheric neutrino deficits are explained through two 
independent sets of parameters and the third generation does not play any more significant 
role than to provide this set. In contrast, in case {ii), the presence of a large mixing between 



^Allowed regions of neutrino oscillation parameters solving the solar neutrino problem in sce- 
nario (i) have been found in ref. [I6|. 
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the first two generations can influence the values of A31 and ^13 allowed by the MSW solution. 
The depletion of the solar neutrino is caused in this case not only by the resonant conversion 
of the neutrinos inside the Sun but also by their oscillations in vacuum arising due to large 
mixing between the first two generations. Thus one expects non-trivial departure compared 
to the case involving only two generations. In this letter we study quantitatively the case (ii). 

The MSW mechanism in the case of three generations is more involved because of the 
presence of the two independent (mass)^ differences and two additional angles as compared 
with two generations. However in a number of situations one can effectively reduce 
the problem to that of two generations. Case (ii) which requires vastly different values for 
IA31I and IA21I falls under this category. One of the mass eigenstates is expected to remain 
unaffected by matter and the other two to undergo a resonant transition. This can be seen 



analytically by means of an approximation [|T9|. Let us assume that mixing between the 
second and the third generation is very small and neglect it altogether. The evolution of 
neutrino states inside the Sun is governed by a (mass)^ matrix which can be conveniently 
expressed in flavor basis as follows: 



Mi 



R 



12l,(7i2;JXl3lt/i3 



V 



^ A 0^ 



R 



R 



T 



13l'^13;-fl'12l'^12 








(1) 




A21 
A31 

where we have subtracted a term m^J proportional to the identity matrix from the actual 
(mass)^ matrix. Rij is rotation in the (ij) plane described by the mixing angle 6ij. A = 
2\/2G pnE with n and E denoting the electron density inside the Sun and the neutrino 
energy respectively. One could exactly diagonalize the upper 2x2 block of this matrix by 
a rotation in the 12 plane by an angle 6*12^- 



{M'^)'' = RU0i2A)MlRu{ei2Ay, 



(2) 



tan^ 



12A 



2(M2 
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We are interested in the limit IA21I ^ l^sil and A. In this case, the 9i2a depends very 
mildly on A: 

Oi2A ~ O12 + e; with e = sin 2^12— —. (4) 

In addition, the matrix (M^)^ also assumes a simple form when non leading terms of 0(e) 
are neglected: 



+ 0(e). (5) 



sfgAgi + Acl2 S13C13A31 

A21 

S13C13A31 A31C13 

where Cij = cosOij; Sij = sin Oij. It is seen from the above equation, that in the limit e going 
to zero, the second generation decouples and MSW resonance could occur between the first 
and the third generation. The resonance condition however differ in this case from the two 
generation case. As follows from eq.(5), the resonance occurs if 

A cos2(ei2) = A31 cos(20i3). (6) 

For the range of A21 that is of interest to us from the point of view of the atmospheric 
neutrino deficit, e < 10^^ and to a good approximation a resonance occurs between only 
two of the neutrinos. The third mass eigenstate remains more or less independent of A. 
This argument is confirmed by the exact calculation of the variation of the effective neutrino 
masses in matter with the solar density displayed in fig. 1 for typical values of the parameters. 

We shall work with eq.(5) with 0(e) terms neglected and comment upon validity of such 
an assumption later on. The probability Pu^u^ fo^' the Uf, produced in the Sun to remain i/g at 
the detector follows immediately using eq.(5) and the standard MSW picture |jTl[. Initially, 



the would be a mixture of three mass eigenstates determined by the form of the mixing 
matrix R at the point of production. One of the mass eigenstates remains independent of A 
while the other two change and undergo resonant transition. Let X denote the probability 
of non-adiabatic transition between the two transforming states. The survival probability is 
then given by 



Pu... = \ (4(1 + (1 - 2X) cos(2^i3) cos(2^^i3Ao)) + 2^2) • (7) 

As follows from eq.(4), the angle 6*12^ is mildly dependent upon A and we have used its 
vacuum value in the above equation. The 6*13^1^ is the value of ^13^1 at the production point. 
This can be determined from eq.(5) 

t«n^ - ^31 Sin(2gi3) 

tan 6^13^0 = -T 77Tr~\ ~a — T- W 

A31 cos(26'i3) - A0CI2 

The exact expression for X depends upon the assumed variation of the density in the vicinity 
of the resonance layer. We shall use the following form which is valid when this variation is 



linear but which provides a good approximation to the actual situation in most cases 

vrAgi sin^ 2^13 



X = exp 



a X 



(9) 



4 Ecos(2^i3 

The above expression follows from a straightforward generalization of the method used in 



case of purely two generation case |TT|. Unlike the resonance condition, eq.(6), X turns out 



to be independent of the large mixing 612 between the first two generations and retains its 



form given in purely two generation case |TT 



We have performed several numerical tests of the accuracy of the analytical expressions 
for the electron-neutrino survival probability given by eqs.(8). In particular, we have checked 
the accuracy of this expression for several values of the parameters Aij and %. This has been 
done by comparing the results of explicit numerical calculation of the survival probability 
with the corresponding values given by eq.(8). One potential source of discrepancy between 
the values to be compared is the nonequivalence of the two approaches. The analytical 
formula presumes an averaging over all vacuum oscillation lengths. On the other hand, the 
numerically obtained survival probability at the surface of the Sun has still to be averaged 
over two oscillation lengths L12 and L31. It is not clear a priori in which order this averaging 
should be performed numerically. We have obtained good agreement between the analytical 
and numerical results when we averaged first over the shorter oscillation length L12 and then 
over the longer one L31. In this case typically the relative error of the analytical formula 
does not exceed several percent. 



Having checked numerically the accuracy of the analytical formula we have proceeded 
with the analysis of the results from the solar neutrino experiments in terms of three neutrino 
oscillations in matter. We have used the latest available data from the four solar neutrino 
experiments taking data at present. For the ratio of the experimentally measured to the 
theoretically predicted event rate in the chlorine experiment |^ we have used the value 

Rci = 0.28 ± 0.03. (10) 

The relevant ratio determined from the published results of the Ga — Ge experiments SAGE 
§ and GALLEX [| is 

i?SAGE = 0.44 ±0.21. (11) 



i?GALLEX = 0.66 ± 0.12. (12) 

We have used also the latest result from the scattering experiment conducted by the 
Kamiokande collaboration ^ 

RKamioka = 0.49 ± 0.08. (13) 



In our calculations we have used the solar model |]20| which takes into account the diffusion of 
hellium and is in excellent agreement with the data from helioseizmology. The restrictions 
on the parameters A31 and ^13 implied by the solar neutrino experiments depend upon 
the value of 612- We have chosen values of 612 allowed by the combined search for the z/^ 
oscillations in atmospheric as well as accelerator neutrino experiments. These values, as 



determined in ref. JT^, fall in the range 0.35 < sin^ 26'i2 < 0.7. They were determined for 



the case of two generations. But the presence of the third generation is not expected to 
significantly influence this determination as long as the u^- — and z/g — mixings are small 
as assumed here. We shall use two representative values of 612 which fall in the allowed 
band in the following. We also do not expect the results of the analysis of the solar neutrino 
deficit with the effect of the Earth on the oscillations of two neutrinos leading to unobserved 
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seasonal and day-night variations |21|] to change considerably with the inclusion of the third 
neutrino. The mixing angle ^13 relevant for the solar neutrino oscillations is too small to 
affect the neutrino oscillations in the Earth because of the much smaller dimensions of the 
Earth as compared with that of the Sun. Also the oscillation length corresponding to A31 
is much larger then the diameter of the Earth. 

The neutrino survival probability for each set of neutrino oscillation parameters has been 
computed by numerically finding the resonant point and substituting in eq.9 the correspond- 
ing value of the logarithmic derivative of the electron number density. With the neutrino 
survival probabilities so obtained we have calculated the corresponding event rates in each 
detector and have compared them with the ones predicted in the standard solar model. The 
allowed regions of parameters for which the neutrino conversion hypothesis under consider- 
ation cannot be rejected at certain c.l. are determined by minimizing the function: 

i=Cl,Ga,H20 \ '^^ / 

They are shown in fig. 2 for two different values of the large mixing angle 612- Note that with 
diminishing the mixing angle Ou the allowed region of parameters A31 - sin^ 26'i3 converges 



to the one determined in the two-neutrino oscillation scenario ||22|. As a less trivial result 
from the above calculation it follows that a relatively large part of the adiabatic region of 
parameters, i.e. the horizontal branch of the MSW triangle between sin^ 26*13 = 0.07 and 
sin^ 26*13 = 0.9 is allowed at the 95 % c.l. for sin^ 26^13 = 0.65. The mixed solution, i.e. 
the upper part of the hypotenuse of the MSW triangle, is also allowed at the 90 % c.l. for 
sin^ 26^12 = 0.65 and at the 95 % c.l. for sin^2^i2 = 0.45. This result is characteristic of 
the three-neutrino mixing scheme we are considering here whereas in the "standard" two- 
neutrino oscillation scenario the adiabatic and mixed solutions are ruled out at the 95 % c.l. 



from the combined analysis of the results of the four solar neutrino experiments ||2^ . 

These results change considerably when we include in our analysis the recoil electron 
energy spectrum data obtained by the Kamiokande-II collaboration. Instead of comparing 
only the mean value of the suppression factor RH2O with the predicted ones for each set of 



neutrino oscillation parameters, we have deleted the corresponding summand from eq.(13) 
and have added instead the function for the 12 recoil electron energy bins as given in ||23|| . 
We end up with the plot shown in fig. 3. In this case only a part of the so called "mixed" 
solution survives but most of the adiabatic region is ruled out at the 95 % c.l.. However, the 
allowed region of parameters is still different from the corresponding one in the two-neutrino 
oscillation case. 

Solar neutrino experiments that are in operation, as well as proposed ones that are being 
developed will make it possible to test the proposed three-neutrino oscillation solution of 
the solar neutrino problem. For sin^ 29i2 = 0.65 and sin^ 2^13 = 5 x 10~^ and A31 = 8 x 10~^ 
eV^ the signals in Ga-Ge detectors should be close to 82 SNU, in chlorine detectors about 
2.52 SNU and the signal in water Cherenkov-detectors should be about 0.33 of the predicted 
one within the standard solar model. For sin^ 2^13 = 6 x 10^'^ and A31 = 5 x 10^^ eV^ the 
signals in Ga-Ge and chlorine detectors should be correspondingly 53 SNU and 2.20 SNU 
and the suppression of the signal in u^e scattering experiments should be 0.45. Finally, for 
sin^ 2^13 = 0.7 and A31 = 10"^ eV^ the signals in these detectors should be 55 SNU, 2.28 SNU 
with a suppression of 0.34 of the signal in neutrino-electron scattering experiments. The pp 
neutrino signals alone should be suppressed by factors of 0.66, 0.63 and 0.44 correspondingly, 
whereas 0.862 MeV ^Be neutrino signals should be suppressed by factors 0.67, 0.47 and 0.44 



respectively. Thus by comparing the signals in Borexino SNO [|25[, ICARUS ||26| and 



Superkamiokande it will be possible to pinpoint the allowed region of parameters. 



The proposed hellium detector of solar neutrinos HELLAZ p8[ is expected to measure 
the fluxes of both pp and ''-Be neutrinos. Moreover, the spectrum of pp neutrinos will possibly 
be measured too. As shown in p9[ this spectrum does not depend on details of the solar 
model and its shape can be predicted with a great precision. Any significant deviation from 
the standard shape will be an unequivocal evidence in favour of neutrino conversion taking 
place either in the Sun, between the Sun and the Earth and/or in the Earth. 

The scenario we have considered requires a hierarchy IA21I ^ I ^31 1- Such hierarchy 
does not occur in seesaw models based on grand unified groups such as 5*0(10) if the right 
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handed neutrino masses are assumed to be flavor independent. However there exist other 
mechanisms for neutrino mass generation which could lead to the required values for the A jj . 
A concrete example is provided by the model of Zee which in fact comes very close to 
the present scenario. This model contains two hierarchical Aij. The larger one describes the 
Uf, — z/^ oscillations while the z/g — Uj. oscillations are of longer wavelength and are controlled 
by the other (mass)^ difference ||3l]|. The model contains only left handed neutrinos whose 
majorana masses are radiatively generated and are described by the following mass matrix 



= mo 



a cos a 
a sin a 
cos a sin a 



(15) 



where mo, c and a are parameters deflned by Wolfenstein [^]. In the limit a going to zero, 
one has a massless state and two with masses ±mo. We need to identify the massless state 



with "z/^". To flrst order in a, one has [0 



mi = mo(l + l/2o"sin2a) 

m2 = —cr sin 2a 

ms = — mo(l — l/2(Tsin2a) 



(16) 
(17) 
(18) 



Hence, A21 ~ mg ; A31 ^ 2mocrsin2Q;. Thus our scenario gets realized by choosing mg ~ 
10~3 (eV)2 and 2a sin 2a ^ -IQ-^. 

The mixing pattern predicted in the model also comes close to the one required here. 
The mixing matrix is given in the limit a going to zero by: 



UMyU'^ = diag.{mi,m2,rn3] 



(19) 



1 / -\/2 cos a sin a 1/ cos a 
sin a — cos a 1 / -\/2 sin a 
1/^2 -1/^2 

This coincides with the mixing matrix that we have used namely. 



(20) 



U' = i?12(^12)i?13(^13), (21) 

if one identifies 612 = tc — a and ^13 = 7? + n/A. a and fience 612 is not fixed within the 
model but would be required to lie in the range appropriate for solving the atmospheric 
neutrino problem. In contrast, sin^ ^13 is already fixed around |Q its maximal value. This 
value is seen to be allowed by the restrictions coming from the solar neutrino experiments 
as displayed in figs. (2, 3). Thus, in addition to having the required hierarchy in masses, the 
mixing angles in the model are also consistent with the solar and atmospheric neutrino data. 

Conclusions. We have attempted to describe the observed deficits in the solar and 
atmospheric neutrinos in terms of the neutrino oscillations involving three generations. The 
deficits in the flux of the low energy atmospheric z/^ is not as clearly established |^ as the 
one observed in case of the solar neutrinos. However if such a deficit gets firmly established 
then the three generation scenario considered here would provide an interesting mechanism 
for understanding both deficits in a coherent manner. As discussed above, the conventionally 
employed two generation MSW picture is still approximately applicable in the present case 
but the large mixing 612 changes the allowed region of parameters in a significant manner. 
Moreover, Uf, — oscillations with relatively large (mass)^ difference envisaged here can 
be studied in laboratory as well. Future experiments may confirm or rule out the present 
scenario. 
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Figure Captions 

Fig.l Masses of three flavors of neutrinos as a function of A expressed in units of 10~^ eV^ 
for tlie values of tlie vacuum parameters: ml = 0; m2 = 500; m| = 1 in units of 10"^ eV^; 
Oi2 = 40°; = 5° and ^23 = 0. 

Fig. 2 Allowed region of parameters sin^(2^i3) and Amsi at the 90 % c.l. (solid line) and 95 % 
c.l. (dashed line) from the analysis of the results of the four solar neutrino experiments. The 
values of the mixing angle 612 have been chosen from the region allowed by the analysis of 
the deficit of atmospheric in terms of oscillations of two neutrino flavours in the Earth , 
sin2(2^i2) = 0.65 (fig.2a) and sin2(2^i2) = 0.45 (fig. 2b). 

Fig. 3 The same as in fig. 2 with the recoil electron energy spectrum as measured by the 
Kamiokande-II collaboration taken into account. 
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